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Abstract — This paper presents a combined Monte Carlo and 
analytic approach to the calculation of the pixel-to-pixel 
distribution of proton-induced damage in a HgCdTe sensor array 
and compares the results to measured dark current distributions 
after damage by 63 MeV protons. The moments of the 
Coulombic, nuclear elastic and nuclear inelastic damage 
distributions were extracted from Monte Carlo simulations and 
combined to form a damage distribution using the analytic 
techniques first described in [1]. The calculations show that the 
high energy recoils from the nuclear inelastic reactions 
(calculated using the Monte Carlo code MCNPX [2]) produce a 
pronounced skewing of the damage energy distribution. While 
the nuclear elastic component (also calculated using the MCNPX) 
contributes only a small fraction of the total nonionizing damage 
energy, its inclusion in the shape of the damage across the array is 
significant. The Coulombic contribution was calculated using 
MRED [3-5], a Geant4 [4,6] application. The comparison with 
the dark current distribution strongly suggests that mechanisms 
which are not linearly correlated with nonionizing damage 
produced according to collision kinematics are responsible for the 
observed dark current increases. This has important implications 
for the process of predicting the on-orbit dark current response of 
the HgCdTe sensor array. 

I. INTRODUCTION 

ANY emerging space astronomy programs will perform 
their science using infrared detectors in order to study 
the early Universe as well as Earth and planetary sciences, and 
the infrared bands are also important in other civil and military 
applications. Although we have observed hot pixel formation 
in proton-irradiated Rockwell IR hybrid detectors to be used in 
the James Webb Space Telescope (JWST) [ 8 ], we do not as 
yet understand the mechanism producing the hot pixels in 
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HgCdTe. As a result we are unable to predict hot pixel 
formation on orbit. The purpose of this paper is to predict the 
proton-induced displacement damage distributions in 
Hgo. 7 Cdo. 3 Te detector arrays based on collision kinematics in 
order to see if they predict the observed dark current 
distribution. 

In the case of Si sensors (including charge couple devices, 
active pixel sensors and charge injection devices) 
measurements show that the dark current distributions are 
often well explained by the damage distributions calculated 
based on collision kinematics [1,9-12]. Damage distributions 
were first calculated analytically by Marshall et al. in 1990 
with good agreement obtained for dark current distributions 
produced by 12 MeV protons in Si charge injection devices. 
At 63 MeV the data indicated less variance in the measured 
distribution than in the damage energy calculation, a result also 
found by Hopkinson et al. at 100 MeV in Si charge coupled 
devices (CCDs). Using the Monte Carlo code CUPID [13], 
Dale et al. showed that this result followed because the recoil 
ranges were comparable to the size of the dark current 
sensitive volume. In the limit of bulk material, both the 
analytic and the CUPID Monte Carlo approaches are in good 
agreement. As sensitive volumes shrink and incident proton 
energies increase, the ranges of the spallation recoil fragments 
approach the smallest dimension of the microvolume, and the 
pixel-to-pixel damage variance are best calculated using 
methods which track the damage deposition along the recoil 
atom pathlengths. In this regime, a Monte Carlo approach is 
well suited to describe the damage energy distribution. 

Nevertheless, in some cases the Si dark current distributions 
cannot be described using collision kinematics alone. This has 
been attributed to hot pixel formation from electric field 
enhanced emission [e.g. 11, 14-17]. It is important to 
distinguish the two scenarios because in order to predict the 
hot pixel populations and assess dark signal nonuniformity, 
one needs to understand whether the underlying mechanisms 
are due to extraordinarily large damage regions from inelastic 
reactions in the pixel or by more ordinary damage in the 
presence of electric field enhanced emission which appears to 
follow from small damage regions in veiy small microvolumes 
filling only a tiny fraction of the pixel’s volume. In the later 
case the hot pixel population may be expected to follow the 
Coulomb cross section [14,17], 



To be published in Transactions on Nuclear Science (TNS) Special Edition, Institute of Electrical and Electronics Engineers 
(IEEE), June 2007. 


1 



? 


II. The Experiment 

We use a previously measured dark current distribution of a 
Rockwell H2RG which is a hybrid device with a 2k x 2k 
format and 18 micron pixel pitch. It incorporates a software 
configurable silicon readout circuit bump bonded to a HgCdTe 
detector array optimized for the JWST NIR (0.6 ~ 5.5 pm) 
spectral range. Details of the experiment can be found in [ 8 ] 
but key details are provided here. 

An engineering grade device was employed and displayed a 
number of cosmetic defects and ! hof pixels that did not meet 
the stringent JWST flight focal plane array (FPA) operability 
requirements going into our test. A subset of -266,000 pixels 
were extracted that were deemed to be “good” pixels as will be 
described in the next paragraph. The detector was held at 37 
K and irradiated with 63 MeV protons to a 5 krad(Si) level 
which corresponds to a fluence of 3.7 x 10 10 cm' 2 . The dewar 
was maintained at temperature while being transported from 
the UC Davis Crocker Nuclear Laboratory back to the NASA 
Ames laboratory, and measurements were taken after residual 
radioactivity from the proton exposure had mostly decayed. 
Residual activity and cosmic ray effects were filtered out as 
described in [ 8 ]. 

As noted in [ 8 ], because of the very low dark current levels 
in these devices, the number of 'good' pixels was difficult to 
precisely quantify in the presence of read noise, cosmic-ray 
events and induced radioactivity in the cryostat. In order to 
isolate and remove engineering defects from the test data, a 
long series of darks was obtained in a clear environment with 
each device prior to the proton exposure. The distribution of 
dark currents estimated from the longest integration for each 
sensor chip assembly (SCA) was examined to determine the 
width of the peak in the histogram and thus the accuracy to 
which we could estimate dark current in these devices. This 
distribution width arises from a combination of system read 
noise and the true non-uniformity in dark current among the 
best pixels in the device. From the observed distributions it 
appeared that read noise was the greatest contributor to this 
width for each test SCA. We chose to set a dark current 
threshold for each long-integration frame in the dark current 
series at the median-pixel dark current plus the distribution 
width of the longest integration frame, as given by its full- 
width-half-maximum. By applying this threshold in each of the 
long-integration frames and removing from our select 
population any pixels that exceed the threshold in any of the 
frames, we derived a subpopulation of pixels that never exceed 
the dark current threshold over many hours of data collection. 
Post-irradiation data reduction is limited to this very 
conservative subpopulation of selected pixels. Note that this 
method is guaranteed to misidentify a number of good pixels 
as bad, since it makes no attempt to correct for the effects of 
cosmic ray events that occur during the long integration series, 
and because read noise will occasionally push marginally good 
pixels above the threshold. We consider this is not a problem 
in this work as we still had a sufficient population of pixels 
remaining to study the dark current distributions. 


III. Damage Energy Distribution Calculations 

In this section we follow the method described in [1] for Si, 
but with modifications, to calculate the damage energy 
distribution for 63 MeV protons on Hgo. 7 Cdo. 3 Te. The first 
step in the calculation of damage energy distributions is to 
calculate the interaction cross sections (a), as well as the mean 
damage energy, (p) and the associated variance, (cr 2 ) of the 
Hgo.7Cdo.3Te damage energy distributions due to Coulombic, 
nuclear elastic and nuclear inelastic interactions respectively. 
These means and variances correspond to the probability 
density function (pdf) governing the likelihood of a particular 
recoil energy resulting from a given proton- Hgo. 7 Cdo. 3 Te 
interaction. 

Underlying assumptions include the fact that interaction 
mechanisms are random in nature, and for our purpose of 
examining the possible correlation of device dark current with 
displacement damage, we consider all non-ionizing energy; as 
with the related Nonionizing Energy Loss (NIEL) energy 
deposition rate. To achieve this, we use the Lindhard partition 
[18] as applied to the calculated recoil energy for a given ion 
and energy. We considered that the nonionizing energy 
associated with a recoil atom was deposited in the pixel in 
which it originated. In the next section we describe how the 
damage energy contributions from Coulomb scattering, nuclear 
elastic scattering, and nuclear inelastic scattering are 
determined, and show that each mechanism has its own 
probability density function (pdf) describing the probability of 
an individual pixel receiving a given amount of damage. As 
the three mechanisms are independent, we then show how 
relatively straight forward statistical tools can combine these 
pdfs in order to describe the damage distribution throughout 
the array; taking into account the pixel geometry, material 
composition, proton energy and proton fluence. Our treatment 
exercises these tools for the case of comparing a measured 
dark current distribution with a predicted damage energy 
distribution to test the hypothesis that there may be a linear 
correlation. 

A. Material and Recoil Spectrum Parameters 

The Hgo.7Cdo.3Te has a density of 7.41 g/cnf and a gram 
molecular weight of 151 g. The JWST H2RG pixel area is 18 
pm by 18 pm and the HgCdTe layer is -10 pm thick. From 
measurements on other HgCdTe detectors we expect the 
diffusion length to be at least 10 pm, so dark current should be 
collected from the entire pixel volume, even though the diode 
junction occupies a small fraction (- 10 %) of the pixel volume. 
In this geometry, the volume and density are large enough so 
that we do not expect the ranges of the recoils to be long 
compared to the volume dimensions so the analytic approach 
should offer a valid approximation. 

The results for the recoil spectrum parameters are shown in 
Table 1, along with the MRED (Monte Carlo Radiative Energy 
Deposition) [3-5] results for the Coulomb parameters in the 
units indicated. MRED is a Geant4 based tool that employs a 
shielded Coulomb potential. In the future, we plan to develop 
MRED for a more comprehensive and general solution but for 
now we are using it for a “point” solution for the case of 3.7 x 
10 10 cm ' 2 63 MeV protons on the pixel geometry cited (i.e. 1 .2 
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x 10 5 protons/pixel). This proton exposure results in 5 
krad(Si). 


Table 1 

Proton Nuclear Reactions on HgCdTe x=0.3 


Proton 

Energy 

(MeV) 

Cross 

Section 

(barns) 

Mean 

Recoil 

Energy 

(MeV) 

Mean 

Damage 

Energy 

(MeV) 

Variance of 
Damage 
Energy 
(MeV) 5 

Nuclear Elastic Reactions 

63 

1.566 

0.0367 

0.0236 

5.56x1 O' 3 

Nuclear Inelastic Reactions 

63 

1.599 

0.365 

0.204 

2.87 xIO" 2 

Coulombic Interactions 

63 



0.734* 

9.60 x 10‘ 13 ** 


*MeV per pixel with incident 63 MeV proton fluence of 3.7 x 10 10 cm' 2 . 
** MeV 2 per pixel with a 63 MeV proton fluence of 3.7 x 10 10 cm' 2 . 


We arrived at the distribution describing Coulomb damage 
based on the mean damage energy deposited in a large number 
of Monte Carlo runs (using MRED [3-5]]) for the pixel 
geometry and proton fluence of 3.7 x 10 10 cm' 2 63 MeV 
protons. Using 1 00 separate simulations of the mean damage, 
we calculated the sample mean presented in the Table. It is 
interesting to note that the mean damage energy per pixel 
calculated using the NIEL value in Fig. 1 (based on the 
analytic Ziegler Biersack Littmark (ZBL) method [20]) agrees 
to within 17% with the value obtained using the MRED Monte 
Carlo runs. 

Calculation of the sample variance for the 100 simulations 
yields the surprisingly small value listed in Table 1. Note that 
each simulation incorporated a new random number seed, and 
each case was run for exactly 1.20 x 10 5 protons/pixel. 
Consideration of our test condition and inclusion of the 
variance associated with the Poisson probability describing the 
incident particle fluence associated with an average of 1.2 x 
10 5 protons/pixel yields a much larger variance of 4.49 x 10” 6 
MeV 2 . The variance used in the following calculations is 
therefore very heavily dominated by this Poisson contribution, 
and for such a large incident fluence we have confidence in 
assuming a Gaussian form for the distribution of Coulomb 
damage throughout the array. Inspection of the variance 
indicates very little difference in the Coulomb damage from 
pixel-to-pixel. This follows from the fact that the mean energy 
imparted per Coulomb scattering event is on the order of the 
displacement threshold (-20 eV), and this results in the 
generation of isolated Frenkel pairs (and possibly point defects 
in HgCdTe). So the Coulomb damage is spread very evenly 
throughout a pixel and hence the pixel-to-pixel uniformity is 
expected. This parallels the treatment of proton damage from 
Coulomb scattering in Si in [19,20 and references therein]. 

The nuclear elastic component (calculated using the Monte 
Carlo code MCNPX [2]) has an almost negligible effect on the 


NIEL (Fig. 1) as it is typically an order of magnitude below 
the dominant term across the entire energy range. But as we 
will demonstrate later, consideration of the nuclear elastic 
contributions has important consequences for the shape of the 
composite damage distribution. We can see from the table that 
the mean recoil energy and variance are an order of magnitude 
below that of the inelastic interactions, but their cross sections 
are essentially equal. The difference in damage energies per 
event gives rise to the lower NIEL relative to the nuclear 
inelastic events, and this is also illustrated in Fig. 2 which 
shows the comparison of damage energy distributions for 
nuclear events generated in the MCNPX calculations. 



Fig. 1 Proton NIEL in Hgo. 7 Cdo. 3 Te as calculated in [21]. Note that the 
NIEL is very insensitive to the exact stoichiometry. 



Fig. 2 Comparison of the damage energy distributions for the nuclear 
elastic and nuclear inelastic interactions at 63 MeV. Note the abscissa is on a 
log scale to represent the two distributions on the same figure. 

The method used to compute the nuclear contribution to 
NIEL and the associated variance is based on the thin target 
approximation using MCNPX and a methodology developed 
by Jun [22]. A thin cylindrical disc of the material of interest 
with a normalized density of 0.01 atoms/bam-cm was 
modeled, and a simulated pencil beam of protons penetrates 
the material. Using the damage energy tally, the history tape 
written by MCNPX was analyzed to calculate the mean 
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damage energy per source particle, T dam , which is the 

nonionzing portion of energy deposited (i.e. after application 
of the Lindhard partition function). Then, NIEL is calculated 
by: 


5 


NIEL ~~ 



T 


dam 


(*.■*) 


( 1 ) 


where N A is Avogadro’s number, A is the gram atomic weight 
of the target material, N v is the atom density and x is the target 
thickness. By using MCNPX, we were able to compute the 
nuclear contributions to the proton NIEL for each material and 
then superimpose them to arrive at the NIEL for the 
Hgo.7Cdo.3Te. The production of displaced atoms is dominated 
by the Coulombic interactions below 10 MeV, while the 
nuclear collisions (particularly the nuclear inelastic) take over 
at energies exceeding ~50 MeV. For the Coulombic NIEL 
shown in Fig. 1, the calculation was done analytically using the 
ZBL method. This analytic approach provides the mean, but as 
previously mentioned the MRED calculations were necessary 
to determine the variance for Coulomb scattering events. Mean 
values and their variances for nuclear elastic and inelastic 
damage distributions were determined using the history tape 
out of MCNPX for the distributions shown in Fig. 2, and later 
we refer to these as the moments of the single event probability 
density functions for describing nuclear collision damage. 

B. Calculated Damage Distributions 

The cross sections for nuclear scattering events to produce 
damage listed in Table 1, combined with the proton fluence 
chosen to represent a relevant JWST mission exposure, result 
in only a few events in the H2RG pixel volume. We define r\ 
as the average number of recoils per pixel and it is given by 

rj = a ■ ■ p • V pjxel ■ N A / A (2) 


nuclear inelastic collisions are present their importance in the 
damage of a given pixel plays a dominant role. 

Since the damage from Coulombic, nuclear elastic and 
nuclear inelastic events are each independent and random 
variables, the statistical description of the cumulative damage 
can be approached by evaluating them independently and then 
combining by convolution. The functional form of the 
Coulomb portion is Gaussian arising from the fact that each 
pixel has a very large number of relatively small collisions. 
Inspection of the functional form of the damage distributions 
of Fig. 2 reveals an asymmetric distribution which is skewed 
towards higher energies (note the abscissa is on a log scale to 
represent the two distributions on the same figure). After 
considering several candidate distributions, we determined that 
the two parameter gamma function is well suited for describing 
both the nuclear elastic and nuclear inelastic damage functions 
(as in [1]). The gamma function is expressed as r(p*A.) where 
jm is the mean recoil energy and X = p/a 2 where a 2 is the 
associated variance. This result affords a very convenient 
approach for the techniques we use to combine distributions in 
that the convolution of two gamma functions yields a new 
gamma function. This fact, along with the result that 
distribution means and variances add linearly under 
convolution, affords a concise approach to this analytic 
technique. In addition, since the Gaussian function is a special 
case of the gamma function, their combination under 
convolution results in a gamma function and this result means 
we have a simplified path to combining all three independent 
variables without requiring FFT-based convolution methods. 

Using the technique just described, we proceed to combine 
Coulombic damage with the single event probability function 
(SEpdf) for elastic nuclear reactions. This step reduces to 
defining a new probability density function based on the 
gamma function (with added means and variances) as shown in 
equation below. The probability (P) of having the damage 
energy (DE) is 


where a is the interaction cross section, <L> i s the proton 
fluence, p is the HgCdTe density, V is the pixel volume, N A is 
Avogadro’s number and A is the gram molecular weight of 
Hgo.7Cdo.3Te. We find that the average number of nuclear 
elastic and inelastic recoils per pixel is 5.7 and 5.8 
respectively, and the maximum number of events expected for 
any pixel is -20 for our pixel population. Also note that only 
-800 pixels out of -266,000 have no nuclear elastic 
interactions, and similarly -800 pixels have no inelastic 
collisions. As these are independent variables, the population 
of pixels having only Coulomb damage is given by the product 
of the Poisson probability of no nuclear elastics and the 
probability of no nuclear inelastics which yields only -2 
pixels! From this we see that even though the total damage 
from nuclear elastic processes in small, the consideration of 
this damage is very important in arriving at the appropriate 
distribution describing pixels with no nuclear inelastic 
collisions, essentially adding to the tail of the Coulombic 
contribution thereby increasing the variance. Obviously where 


P(DE ) - 


X 

(A ■ DEp r 


e (Z-DE) 

ITT 


( 3 ) 


where and X = |i/cr and r = (i 2 /<r. (This function is also 
normalized to unit area.) For convenience, we note that a 
gamma function is defined as 

co 

T(cc)= \x (a - ]) -e- x dx (4) 

0 

where a > 0 and x is the independent variable. For example, 
the damage probability function (Equation 3) applies to pixels 
that have Coulombic damage and only 1 nuclear elastic recoil 
for which jlx is the sum of the Coulombic mean damage energy 
and the mean damage energy for one elastic recoil, and 
likewise, a 2 is the sum of the Coulomb variance and the 
variance for one nuclear elastic recoil. By extension we 
generate a family of gamma distributions accounting for pixels 
with multiple nuclear elastic recoils by examining the N-fold 
convolution of the SEpdf with itself to get the pdf for a pixel 
with N elastic recoils. This is done for N = 2, 3, 4, etc, up to 
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the maximum N=19 events expected to occur for the given 
population of pixels. The resulting total elastic (combined 
Coulomb plus nuclear elastic) pdfs for a pixel with N nuclear 
elastic recoils is shown in Fig. 3. Note that as N increases the 
mean of the distribution increases, it broadens and also tends 
towards a gaussian - all expected based on the Central Limit 
Theorem. 



0.5 1 1.5 2 

Damage Energy (MeV) 

Fig. 3 The total elastic pdfs for a pixel with N = 0,1, 2,... 21 nuclear elastic 
interactions. (The N = 0 (Coulombic interactions only) is a gaussian curve 
with a maximum value of 29.5.) 

For the sensitive volumes (V p i xe i) an ^ fhience (O) of interest 
the variation in the number of nuclear elastic (or inelastic) 
recoils per pixel is described by a discrete Poisson distribution 
in which the probability P is 


P(X,T]) = 


(e-»«7 T ) 

x\ 


(5) 


where x = 0 , 1 , 2 ,... recoils per pixel and r\ is the average 
number determined according to Equation 2. Using this 
relation, we weight each of the 19 gamma distributions with 
the Poisson probability for x (nuclear elastic) recoils in a pixel, 
and then sum them to arrive at the combined Coulomb and 
nuclear elastic (which we will refer to as the combined elastic) 
damage energy distribution across the array. Note that XP(x,t|) 
= 1 . 

In Fig. 4 we present the Poisson weighted pdfs as well as the 
sum in graphical form, and using standard statistical 
techniques we calculate the mean and variance of this resulting 
combined elastic function to be 0.87 MeV and 0.19 MeV 2 . 
Note that the Poisson weighting and superposition results in a 
functional form which is not described exactly by a gamma 
function, but due to the ease of working with gamma 
distributions for the next step of the analysis to incorporate 
nuclear inelastic damage, we approximate the resulting 
combined elastic distribution with a gamma function having 
the same mean and variance. This function is also shown in 
Fig. 4. The comparison shows a reasonably good fit, and we 
consider the tradeoffs associated with this assumption to favor 
this approximation. 


Note that the variance in the combined elastic damage is 
four orders of magnitude greater than the Poisson dominated 
Coulomb only case. Even though nuclear elastic collisions 
impart only a small fraction of the total NIEL, their inclusion 
in the population of pixels in the array (-800 total) which have 
no inelastic damage is necessary. 

Fig. 1 and Table 1 reveal the important implications of the 
nuclear inelastic damage in that it accounts for >60 % of the 
total damage at 63 MeV, and with an average of 5.8 events per 
pixel the variation from pixel-to-pixel is obviously the most 
significant aspect of the distribution of damage throughout the 
array. Using the techniques already described, we combine the 
N-fold convolution of up to 20 inelastic events per pixel with 
the gamma function describing the combined elastic damage 
and apply Poisson weighting per Equation 4 and finally 
perform the summation of all the Poisson weighted 
distributions. 



0 0.5 1 1.5 2 

Damage Energy (IVIeV) 

Fig. 4 Comparison of the combined elastic damage energy 
distribution and a gamma function with the same mean and variance. 
The distributions underneath the two curves are the Poisson weighted 
gamma distributions from Fig. 3 for pixels with N nuclear elastic 
recoils which are summed to arrive at the total elastic damage energy 
distribution. For clarity only N = 0-6, 8,10 are shown. Note the 
reduced variance for the N - 0 (Coulomb only) case. Recall that the 
average number of nuclear elastics is 5.7. 

The calculated damage distribution is shown in Fig. 5, 
which shows the distribution of pixels with elastic damage plus 
0, 1,2, etc. inelastic interactions underneath the sum of each 
contribution. The average number of inelastic collisions per 
pixel was -5, and the maximum number of inelastic recoils 
expected per pixels is 20 for our fluence and pixel population. 
As one can see, skewness in the distribution is primarily 
attributed to the Poisson weighting of the inelastic interactions. 
This function describes the distribution of damage 
corresponding to our proton-induced dark current distribution. 
Note that the shape assumed for the SEpdfs (gamma functions 
in our case) is not very critical. 

It is interesting to compare the damage distribution 
calculated for a high atomic number compound such as 
HgCdTe with the previous Si result [1]. In HgCdTe the effect 
of the inelastics is much more prominent - a higher mean 
energy and a broader peak with a more pronounced tail is 
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observed because the inelastic cross section is 3 times larger, 
the mean recoil damage energy is 30% higher and the variance 
is an order of magnitude larger in HgCdTe as compared to Si. 



Damage Energy (MeV) 


Fig. 5 Calculated damage energy distribution for 266,000 
Hgo.7Cdo.3Te pixels irradiated with 3.7 x 10 10 cm ‘ 2 63 MeV protons. 
Below the summed damage distribution we show the Poisson 
weighted damage distributions for the pixels with only elastic events 
(left most peak barely visible in figure), and the pixels with elastic 
events plus 1 (left peak with maximum at ~1 MeV), 2, 3... 18 inelastic 
interactions per pixel. (Note that a pixel may contain up to 18 
inelastic interactions though only distributions with up to 10 are 
resolved in the plot) 



0.2 0.4 

Dark Current (e-/s) 


Fig. 6 Measured dark current distribution for 266,000 Hgo. 7 Cdo. 3 Te pixels 
before irradiation and after irradiation with 3.7 x 10 10 cm " 2 63 MeV protons. 
The data represent measured absolute values and the negative numbers reflect 
measurement noise in the system. The mean dark current is extremely small. 


IV. Comparison of Damage Distribution with Dark 
Current Histogram 

The measured dark current histogram for 266,000 selected 
pixels is shown in Fig. 6 along with the pre-irradiation 
histogram. In [8], the slight median shift was not investigated, 
and hence was not presented as necessarily real. However, we 


have re-analyzed this data and do find a small shift in the 
median dark current after irradiation to 3.7 x 10 10 cm" 2 63 
MeV protons as indicated in the figure. It is apparent in Fig. 7 
that the calculated damage distribution does not predict the 
measured dark current distribution which indicates that some 
other mechanism than collision kinematics is also responsible 
for the high dark current pixels seen in the data. This makes 
an on-orbit prediction of the dark current based on proton 
results at a single energy (63 MeV in our case) problematic 
since the nonionizing energy loss rate (NIEL) correlation does 
not appear to hold. Measurements at 8 MeV are planned to 
see if the high dark current pixels correlate with the Coulombic 
portion of the NIEL which would be the case if electric field 
enhanced emission via trap-assisted tunneling is responsible as 
has been seen in Si. (Trap-assisted tunneling is known to be 
important in HgCdTe sensors at low temperatures.) If this is 
also the case with the present HgCdTe array then first order 
estimates of hot pixel rates expected on-orbit which are based 
on dark current histograms for the 63 MeV equivalent fluence 
for a given mission will underestimate the number of hot 
pixels. 
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Fig. 7 Comparison of measured dark current histogram with 
damage energy calculation. The high dark current pixels are not the 
result of collision kinematics. 
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V. Discussion and Summary 

In comparison to earlier analytic damage energy distribution 
predictions (e.g. [1]) we find that the role of the inelastics in 
producing the highly damaged pixels is even more dominant 
for the HgCdTe case as compared to Si. The cross section for 
inelastic interactions is larger in the (higher atomic number 
system) HgCdTe by a factor of three and the variance is an 
order of magnitude larger in HgCdTe. 

We have developed a new combined Monte Carlo and 
analytic approach to the calculation of the pixel-to-pixel 
distribution of proton-induced damage in a HgCdTe sensor 
array and compared the results to measured dark current 
distributions after damage by 63 MeV protons. The moments 
of the Coulombic, nuclear elastic and nuclear inelastic damage 
distributions were extracted from Monte Carlo simulations and 
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used to generate single event probability distributions 
describing each class of proton-HgCdTe interaction. These 
were combined as discussed in the text to form a damage 
energy distribution. The calculations show that the high 
energy recoils from the nuclear inelastic reactions (calculated 
using the Monte Carlo code MCNPX [2]) produce a 
pronounced skewing of the damage energy distribution. While 
the nuclear elastic component (also calculated using the 
MCNPX) contributes only a small fraction of the total 
nonionizing damage energy, its inclusion in the shape of the 
damage across the array is significant because its variance 
dominates in those pixels with no inelastic reactions. The 
Coulombic contribution was calculated using MRED [3-5], a 
Geant4 [4,6] application. We noted that only 2 pixels in the 
entire array have only Coulombic interactions. The 
comparison with the 63 MeV proton-induced dark current 
distribution strongly suggests that mechanisms which are not 
linearly correlated with nonionizing damage produced 
according to collision kinematics are responsible for the high 
dark current pixels. Measurements at 8 MeV are planned to 
see if the high dark current pixels correspond to the Coulombic 
portion of the NIEL which would be the case if electric field 
enhanced emission is responsible as has been seen in Si. The 
technique for describing the damage energy distributions is 
extendable to other proton energies, material systems, and 
pixel geometries. 
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